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ABSTRACT 


A  modified  point  mass  mathematical  model  which  incorp>rates  an 
estimate  of  the  yaw  of  repose,  has  been  developed  to  represent  the 
flight  of  a  spin  stabilized,  dynamically  stable,  artillery  shell.  This 
improved  mathematical  model  has  the  desirable  feature  of  representing 
the  effects  of  the  significant  variables  of  yaw  of  repose  and  axial  spin 
along  the  trajectory. 
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TABijE  of  symbols 


Definition  Units 

2 

Axial  moment  of  inertia  lb-ft 

Azimuth  of  line  of  fire  (clockwise  from  north)  mil 

2 

Transverse  moment  of  inertia  lb-ft 

Ballistic  coefficient  for  standard  mass  lb/in^ 

Reference  diameter  of  projectile  ft 

Position  of  projectile  with  respect  to  spherical  ft 


Earth  surface 

Acceleration  due  to  gravity 
Acceleration  due  to  gravity  (surface) 
Total  angular  momentum 
Unit  vector  in  the  direction  of  v 
Spin  damping  moment  coefficient 
Drag  force  coefficient 


ft/sec^ 

ft/sec^ 

lb-ft2-rad/sec 


Yaw  drag  "'rce  coefficient 
Magnus  force  coefficient 
Damping  moment  coefficient 
Lift  force  coefficient 
Overturning  moment  coefficient 
Pitching  force  coefficient 
Magnus  moment  coefficient 
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Term 

1 
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N 

Q 

r 
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t 


w 


X 


A 


Definition 

Lift  factor 
Latitude  of  launch 
Mach  number 
Projectile  mass 

Standard  projectile  mass 

Axial  spin 

Yaw  drag  factor 

Distance  between  center  of  Earth  and  projectile 

Effective  radius  of  Earth 

Time 

Velocity  of  projectile  with  respect  to  ground 

Velocity  of  the  projectile  with  respect  to  air 

Velocity  of  the  air  with  respect  to  ground 

Unit  vector  along  the  longitudinal  axis  of 
the  projectile 

Position  of  the  projectile  with  respect  to  a 
ground-fixed  coordinate  system 

Angle  of  yaw  of  projectile 

Approximation  for  yaw  of  repose 

Coriolis  acceleration  due  to  rotation  of  the 
Earth 


Air  density  as  a  function  of  (where  is 


component  of  E) 


Units 

deg 

lb 

lb 

rad/ sec 

ft 

ft 

sec 

ft/  sec 

ft/sec 

ft/sec 

ft 

rad 

rad 

ft/  sec^ 
lb/ft3 
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T  erm 

a 


Definition 

Angular  velocity  of  the  Farth 

First  derivative  with  respect  to  time 

Second  derivative  with  respect  to  time 


Units 


rad/ sec 

- /  sec 

- /sec 
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«WIC0S  »4ds  WAS  BLUE  THCffiTOR*  WAS  NOT  FltKED 


INTRODUCTION 

The  mathematical  model  for  rigid  body  trajectory  simulation,  as 
reported  in  BRL  Report  No.  1244,  closely  matches  the  results  of 
physical  experiments  of  spin-stabilized  projectiles  over  a  spectrum 
of  test  conditions.  However,  its  usefulness,  as  in  most  rigid  body 
simulations,  is  hindered  by  the  time  required  to  numerically  solve 
the  system  of  differential  equations.  This  report  describes  the  der¬ 
ivation  of  a  mathematical  model  which  will  not  be  as  time  consuming 
to  solve  as  the  rigid  body  system,  but  which  will  represent  as  closely 
as  possible  the  center  of  gravity  motion  of  the  projectile  by  utilizing 
a  force  system,  axial  spin  and  an  estimate  of  the  yaw  of  repose. 


BASIC  LAWS  OF  MOTION  OF  A  RIGID  BODY 


The  frame  of  reference  for  all  vectors  to  be  presented  is  aground- 
fixed,  right-handed  Cartesian  coordinate  system  with  unit  vectors 
(1,  2,  3) .  Assume  the  2  axis  to  be  parallel  to  the  vector  g  ,  and  g 
to  have  the  same  direction  as  X 

Asstime  that  the  body  can  be  considered  a  solid  of  revolution.  Then 
spin-stabilized  projectiles  can  be  oriented  by  choosing  a  unit  vector  x 
along  the  axis  of  symmetry,  pointing  from  tail  to  nose.  N  is  the  mag¬ 
nitude  of  angular  velocity  parallel  to  x.  N  is  positive  if  it  results  from 
a  rotation  causing  a  right-handed  screw  to  advance  in  the  direction  of 
x . 

The  following  set  of  simultaneous  differential  equations  of  motion 
for  a  spin-stabilized  projectile  is  developed  in  BRL  Report  1244. 

The  equation  of  motion  of  the  center  of  mass  is: 

(l.D  u*  -pdZ  (Kp  +  a2KD  )  v  v+£d2KL  [  v  X  (x  X  v)  ] 
m  J-,o  a  m 

-pd^  K„  v  x  +  pd^K  N  {xXv)+g+A 

m  &  “*■  _ F  — 

m 

The  total  angular  momentum  of  the  body  can  be  expressed  as  the 
sum  of  two  vectors  in  the  ground-fixed  coordinate  system: 

(a)  The  angular  momentum  about  x . 

(b)  The  total  angular  momentum  about  an  axis  perpendicular  to  x. 

The  angular  momentum  about  x  can  be  represented  by  the  vector 
A  N  x  and  the  angular  momentum  about  an  axis  perpendicular  to  x  by 
the  vector  B  {  x  X  x ). 
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Let  H  denote  the  total  angular  momentum  of  the  body.  The  vector 
representation  of  H  is: 

(1.2)  H  =  ANx  +  B  (x  X  x) 

The  vector  rate  of  change  of  angular  mo:m,.itum  is  the  sum  of  the 
applied  moments. 

(1.  3)  H  =  ANx  +  ANx  +  B  (x  X  x)  = 

pd^K.rv  (v  X  x)  -  pd^K_Tv  (x  X  x) 

M  -*•  -*■  H  -*•  -*■ 

+  pd^  K^,  N  [  x  X  (x  X  v)  ]  -  pd^  N  v  x 

ESTIMATE  FOR  YAW  OF  REPOSE 

An  estimate  for  the  yaw  of  repose  will  be  derived  from  the  rigid 
body  system  of  differential  equations  of  motion  for  the  purpose  of  rep¬ 
resenting  the  effects  of  yaw. 

The  following  conditions  are  assumed: 

(1)  The  projectile  can  be  represented  sufficiently  well  as  a 
body  of  revolution. 

(2)  The  projectile  i?  dynamically  stable. 

(3)  Initial  yaw  is  assumed  small,  i.  e.  ,  it  has  no  sigi.ificant 
effect  on  the  trajectory. 

Examination  of  equation  (1.  1)  shows  the  magnitude  |  x  X  v  | 
as  being  present  in  the  lift  term  and  the  Magnus  term,  the  lift  term 
being  more  significant.  Examination  of  equation  (1.3)  shows  its 
presence  also  in  the  corresponding  terms  of  H. 

If  I  =  v  /v,  then  |  x  X  I  |  =  sin  a.  To  a  first  order  approximation. 
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this  is  a  ;  however,  for  computational  purposes  it  is  a  more  desirable 
quantity  than  a  .  To  get  the  proper  orientation  of  this  quantity,  the 
following  vector  will  be  defined: 

(2.  1)  a  =  I  X  (x  X  I)  =  x  -  cos  o  I 

->  C  — ►  — ►  — ► 

Obviously,  a_e  .  0.  represents  vector  yaw  *  directed  fromj^ 

toward  x  .  The  effect  of  a  e  on  the  trajectory  is  generally  small  under 
the  assumptions  stated.  Furthermore,  it  will  be  assumed  that  <^e  is 
negligible.  This  implies  a  is  small;  moreover,  the  following  approx¬ 
imations  are  warranted: 

(2.2)  x  =  cos  a  _I 

(2.  3)  x  =  cos  a  J[ 

The  separation  of  H  into  ocmponents  parallel  and  perpendicular  to 
x  yields: 

(2.4)  AN  =  -  pd4K  aNv 

(2.5)  ANx  +  B  (x  X  x)  =  pd^K^vlvX  x) 

-  pd4K  v  (x  X  x)  +  pd^K  N  [x  X  (x  X  v)] 

To  determine  <^e  ,  first  replace  x  and  its  derivatives  in  (1.  1)  and 
(2.  5)  using  equations  (2.  1),  (2.  2)  and  (2.  3). 


*  NOTE:  |f^e|  °>  However,  as  mentioned  earlier,  in  magnitude 

|  a^e  |  is  a  first-order  approximation  for  yaw.  is  in  the  plane 
determined  by  x  andj^  ,  in  the  direction  from^  toward  x.  Hence,  a^e 
will  be  referred  to  as  vector  yaw. 
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The  resulting  equations  are: 

-pd^(K  +  c2K  )v2  1  pd^K  v2  a.  pd^K„v  cos  a  I 

JJ  *-J  .  Li  O 

(2.  6)  b  =  _ 2- _ S _ +  _ - _ 

m  mm 

3 

+  pd  KfNv  (c_e  x^_)  4  g  +  A 
m 

"  **  2 
(2.7)  AN  cos  o  H  B  cos  a  (£g  x!^)  +  B  cos  a(I  ) 

3  7  4 

=  pd  K^v  (^  Xc^e)  +  K  Nv  [  cos  a  (a_e  +  cos  a_I)-^I] 

4 

-pd  K^v  cos  a  [(a^e+  cos  C_I)X^_  ] 


Cross  multiplication  by  I  of  both  sides  of  equations  (2.  6)  and  (2.  7) 
and,  with  A  negligible  in  comparison  to  g,  the  simultaneous  solution 
of  the  resulting  equations  gives  the  following  for  a  , 

(2.  8)  a  =  { mpd^  K  Nv  cos  a  [I  X  (  u  -  g  +  pd^K  v  cos  a  I  )] 

m 

+  pd  v2  [-  AN  cos  a  (  I^X^)  +  B  cos2  a  [J^X  (_I  X  JI)  ] 

+  pd^K^  v  cos  “  c  _I  ]  }  ^  { p2  d^  N2  v2  cos  a 

+  pd2  Kt  v2  [  pd2  K  ,  v2  +  B  cos  cl  .  I  ]  } 
r  L  M  -*  -*  1  J 
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For  further  substitution  into  (2.8)  the  following  are  required: 

(2. 9)  _I  =  v  /  v 

(2.10) 1=  [v-(v.I)_l]  /  v  =  [  _ix  (  v  XJQ  ]  /v 

(2.11)  I_  =  -  2  (v.  I)v  /v2+  ^(v.  I.)2I./v2 

+  V  /  v  -  (  v  .  v  )^  /  v2  —  (v_  jy  1^/  v 
(2.  12)  v  /  v  -  (**  .  !_)_!./  vs  [_I_X  (  v  X_I.)]  /  v 


(2.  12)  will  considered  negligible  since  v  can  be  approximated 
by  u  ,  and  u  is  essentially  parallel  to  I . 

With  this  assumption, 

(2.  13)  c  =  {-  AK  N  cos  a  (v  X  v)  +  md2  K  N  cos  a  [v  X(u  -  g 
Xj  “*■  X 

+  pd3  K  cos  a  v  /  m)  ]  -K[v~(v.I)l] 

O  J-I  -+■  -+  -> 

2 

[  2  B  cos  a  (v  . 

-V 

fpd3KT  K  v4  +  pd5K_K_N2  cos  a  v2 

iVl  X1  1 

+  B  cos  a  [  (v  .  I^2  -  (v_ .  v)  ]  } 

4 

In  the  denominator  of  equation  (2. 13)  the  v  term  predominates  ; 

therefore,  K  B  cos  a  [  (v  .  I)2  -  (v  .  v)  ]  will  be  considered  negligible. 

The  ratio  of  this  term  to  the  remaining  terms  in  the  denominator  is  of 
the  order  of  (g2/ v^).  In  the  numerator,  [v  -  (v  .  _I)  _l]  is  the  component 
of  v  perpendicular  to  the  projectile  flight  in  the  air  coordinate  system. 
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For  most  spin  stabilized  trajectories  |v  -  (  v  •  I)  I  J  is  no  more  than 
the  magnitude  of  g;  hence,  consider 

2  4  2  2 

(2.  14)  K  [  v  -  {  v  .  I)  I  ]  [  2  B  cos  a  (  v  .  I)  -  pd  K  cos  a  v  ] 

ri 

2  4  2  2 

K.  g  [  2  B  cos  a  (  v  .  I )  -  pd  K  cos  a  v  ] 

For  artillery  shells,  this  term  is  generally  bounded  by 
(2.  15)  KL  |  g  |  2  B  •  s  (lOg)  =  20Kl  BgQ  2 

Now 

(2.  16)  -  AKt  N  {  v  X  u)  +  md2  K  N  [  v  X  (u  -  g)  ] 

2  2  2 

Nv  {-  AK  +  md  K  }  KT  pd  v  +  AK  g  sin  <p 
L.  L  Li  Lj  o 

,  ,  g  •  V 

where  cos  <f>  =  ^ 

g  v 

For  artillery  shells,  it  is  usually  true  that 

(2.17)  {-  AK  +  md2  K  }  K..  Pd2  v2  <<  AK,  g  sin  p 

L  T  L,  ’EY  Li  o 
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t 


Comparison  of  (2.  15}  and  (2.  17)  shows  that  the  ratio  of  (2.  14)  to 
(2.  16)  is  less  than  the  order  of  (20gQ  /Nv).  For  high  spin,  (2.  14)  will 
be  considered  negligible. 


Let 

(2.18)  v  =  u-  w 

(2.19)  v  a  u  -  w 


If  the  effects  of  the  pitching  force  (Kg  term)  and  w  are  assumed 
insignificant,  and  cos  a  can  be  approximated  by  1,  equation  (2.  13)  is 
reduced  to 


S-e  = 


-AK  N  (  v  X  u)  +  md^K  N  [vX(u-g)] 

I  i  -*  "♦  1  -+• 


l”i3KLKMv4  +  ,>d5KFKTN2v2 


(2.20 


UTILIZATION 


The  primary  goal  of  the  development  of  a  e  was  the  acquisition  of  a 
mathematical  model  which  would  incorporate  the  effects  of  yaw,  but 
would  not  require  the  computing  time  of  a  complete  rigid  body  simula¬ 
tion.  The  following  representation  was  devised  to  incorporate  <^e  in  a 
modified  point-mass  mathematical  model.  This  representation  includes 
auxiliary  equations  necessary  for  the  numerical  solution  of  the  differen¬ 
tial  equation  of  motion  of  the  center  of  mass. 

The  equation  of  motion  of  the  center  of  mass  is: 


pm 
r  s 

{K  +  K  [  Qa 

144  C  m 

u  jj  e 

o  a 

s 

+  P<*2 

Kl  V2  Ue  +  g  +  A 

m 

KFNQ(“ex  Z) 

m 


where:  C  =  ballistic  coefficient  for  standard  mass 
s 

1  =  lift  factor 
m  =  projectile  mass 
m^  =  standard  projectile  mass 

Q  =  yaw  drag  factor 


*Note:  If  Q  and  1  are  set  to  zero  (0)  in  equation  (3.  1),  this  system 
reduces  to  the  classical  point  mass  equations  of  motion. 
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The  axial  spin  is: 


(3.2)  N 


»t  ,4 
pd 


Ka  Nv 
A 


The  approximation  for  the  yaw  of  repose  is: 

<3.31  2  _  =  (  a  -  a  )  (v  X  u)  -  a  (  v  X  g) 
-*e  b  a 


where:  a  = 
a 


akt  n 

i-J 

pd3KLKMv4  +  ('■’S-Kt*2'’2 


m  KtN 

pd  klkmv4+  p^kfktn2v2 


The  velocity  of  the  projectile  with  respect  to  air  is: 

(3.4)  v  =  u  -  w 

The  position  of  the  projectile  with  respect  to  ground-fixed  Cartesian 
coordinate  system  is: 


(3.  5)  X  =  C  u  dt 

Jo 

The  approximation  for  the  position  of  the  projectile  with  respect  to 
spherical  Earth  surface  is: 


(3.6)  E  = 


X2  +  R  -  <R2  -  Xx2)  1/,Z 


where  R  -  effective  radius  of  Earth. 


The  approximation  of  the  force  of  gravity  is: 

fXl  1 

2  1 

(3.7)  g  =  X2  +  R 


r 


r  =[  XL2  +  (X2  +  R)2  +X32]  1/2 

where  =  value  of  gravity  at  point  of  launch 

r  =  distance  between  center  of  Earth  and  projectile 
The  Coriolis  acceleration  due  to  rotation  of  the  Earth  is  : 


X,  u,  - 

\„u„ 

1  2 

2  3 

X,  u,  + 

Xu, 

1  1 

3  3 

X,u  - 

Xu, 

2  1 

3  2 

For  the  northern  hemisphere,  the  X's  are  defined  by  the  following 
equations.  [  For  the  southern  hemisphere  replace  Lby  (-  L).  ] 

\^=  2  S3  cos  L  sin  AZ 

2  £2  sin  L 

2  £2  cos  L  cos  AZ 

where:  £2  =  Angular  velocity  of  the  Earth  (radians/sec) 

L  =  Latitude  of  launch  point 

AZ=  Azimuth  of  fire  measured  clockwise  from  north 
The  orientation  of  yaw  {  St)  is  the  angle  between  the  plane  containing 
both  v  and  and  a  vertical  plane  containing  v  .  It  is  measured  clock¬ 
wise  from  the  vertical  plane.  If  desired,  T  is  given  by  the  expression: 


(3.9)  tstan 


( v  a  -  v  a  )  v 

_ 1  S3  3  C1 _ 

{V1%-  V2aei)  V  (v2ae3~ 
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The  dimensionless  aerodynamic  coefficients  are  functions  of  many 
dimensionless  power  products,  including  the  dimensionless  shape 
parameters,  Reynolds  number  and  Mach  number.  Aerodynamic  coeffi¬ 
cients  are  defined  with  reference  to  a  specific  set  of  shape  parameters 
and  may  be  expressed  as  functions  of  Mach  number. 


CONCLUSIONS 

A  modified,  point-mass  mathematical  model  has  been  developed 
which  incorporates  an  estimate  of  the  yaw  of  repose.  This  improved 
mathematical  model  has  the  desirable  feature  of  representing  the  effects 
of  the  significant  variables  of  yaw  of  repose  and  axial  spin  along  the  tra¬ 
jectory.  By  the  incorporation  of  this  yaw  of  repose,  the  factors  (ballis¬ 
tic  coefficient,  yaw  drag  factor  and  lift  factor)  used  to  match  empirical 
results  have  been  found  to  vary  little  from  theoretically  determined 
values  over  the  spectrum  of  conditions. 

In  comparisons  between  the  rigid  body  mathematical  model,  the 
most  complete  representation  available,  and  point-mass  representa¬ 
tion,  the  modified  point-mass  model  accounted  for  better  than  90 
percent  of  the  discrepancies  in  the  time -dependent  variables  range, 
height  and  deflection.  These  comparisons  were  made  for  spin- 
stabilized  artillery  rounds  currently  being  used  by  the  U.  S.  Army. 

The  modified  point-mass  model  solution  required  only  approximately 
twice  the  computation  time  of  the  point-mass  solution  while  the  rigid 
body  solution  requires  about  one  hundred  to  one  thousand  times  that 
of  the  point-mass  solution. 
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